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Abstract 

We present the azimuthal correlations in the decay products of the top quark pair created in the 
process gg — >■ tt — )■ bbW + W~ calculated in leading-order QCD. We then discuss deviations in this 
distribution that could arise from New Physics contributions, focusing, in particular, on the case 
of a new heavy scalar coupling to the top quark and produced via the gluon fusion mechanism. 
Within this scenario we consider two cases: one in which the top quark is the only particle with 
a significant coupling to the scalar and another in which there are additional heavy fermions that 
modify the effective coupling of the scalar to gluons. We discuss the effects on the overall cross 
section and azimuthal correlations of the top quark decay products, with particular emphasis on 
the effect of interference between the QCD and New Physics contributions and on the existence of 
CP-nonconserving couplings. 
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I. INTRODUCTION 



As the heaviest known Standard Model (SM) fermion and one which can now be copi- 
ously produced at the CERN LHC, the top quark is an interesting tool in the search for New 
Physics (NP) beyond the SM. That the top quark mass m t is so similar to the vacuum ex- 
pectation value of the SM Higgs boson, (H) = (v/y/2) « 174 GeV, is suggestive of a possible 
relationship with any New Physics at or near the scale of electroweak symmetry breaking 
that is now being probed. Indeed, there has been much theoretical and phenomenological ef- 



fort recently due to evidence seen at the Fermilab Tevatron of a possible "forward-backward 

Li_3|) with models invoking such possibilities as new 
7], extended Higgs sectors new fermions js], and 



asymmetry" in top pair production 
3-5], axi-gluons jg, 



W and Z' bosons 
other exotic scenarios 10]. 

In this paper, we will be concerned not with the angular distributions of the top quarks 
themselves, but with those of their decay products. We put particular emphasis on the 
ability to use these correlations to probe the CP nature of a new scalar resonance, most 
notably in the case that it is not a, CP eigenstate and therefore introduces a new source of 
CP violation. There have been many previous studies looking into the effect of resonances 
of varying spin and CP properties on top pair production, including the impact on the 
invariant mass spectrum and the behavior of various observables sensitive, in particular, 
to CP llM23|. Here, we are interested in how the angular distributions of the complete 



results, including the QCD prediction and its interference with the resonant signal, could be 
detected as a deviation from the prediction of QCD alone. We discuss the results as they 
would appear at the LHC running at a design center-of-mass energy of \fi = 14 TeV and 
note the difficulties that arise in both detection of the signal generally and the interpretation 
of the angular correlations in the CP-violating case. 

The paper is organized as follows: In Sec. [Til we gi ye the results for the azimuthal cor- 
relations predicted by leading-order (LO) QCD. In Sec. IHIl we describe the models used in 
our discussion of New Physics contributions. The effects on the top pair invariant mass dis- 
tribution are then given and discussed in Sec. II VI and the those of the azimuthal correlations 
are given in Sec. [V] We offer concluding remarks in Sec. IVII 



2 




(a) 




(b) 



FIG. 1: Diagrams contributing to the process gg 



w+ 




tt 



(c) 

bbW + W~ at LO in QCD. 



II. QCD PREDICTIONS 

The process gg —> tt — > bbW + W~ is given, to leading-order (LO) in QCD, by the three 
diagrams 1 listed in Figs. [U(a)-[l](c). While bbW + W~ are not the true final-state particles 
- the W bosons will decay and the b quarks will hadronize - we assume that the events 
can be sufficiently reconstructed to be described at this level. In practice, this involves 
significant experimental challenges which vary with the decay modes of the W bosons: in 
the fully hadronic channel, one must worry about the large SM background to the overall 
process bbjjjj in addition to the complicated combinatorics involved with the large number 
of jets, while in the leptonic modes there is difficulty in precise event reconstruction due 
to the missing energy carried away by the neutrinos and the unknown components of their 
momentum along the beam axis (see, e.g., the discussions by the ATLAS and CMS collab- 



orations in Refs. 



24 



251]). There are ways in which these issues can be partially alleviated 



using various constraints, such as requiring the reconstruction of the top quark and VF-boson 



1 Because of the triple gluon vertex in Fig. (Ha) there may in principle be an additional diagram with the 
two external gluons replaced by ghosts, depending of the choice of gauge. We choose a gauge such that 
this diagram is absent but in which the sum over gluon polarizations is modified. This choice is used 
consistently across all calculations throughout the paper. 
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FIG. 2: This defines the geometry relevant to the measurement of the angle <3?. In the it rest frame 
we identify the t quark 3-momentum as the z! direction (which is not to be confused with the beam 
direction in the laboratory frame, often denoted as z). The angle $ is then the relative azimuthal 
angle between the decay planes of the t and t about this axis as indicated. This angle can also 
be understood equivalently as the difference in the azimuthal angles of either the W bosons or b 
quarks about the z' axis. 

masses [20I ] . and we note that our results are independent of the particular W^-boson decay 
modes and can therefore be applied to any reconstructed experimental sample of bbW + W~. 
In what follows, we will assume maximum event reconstruction and selection efficiency for 
illustrative purposes while noting the above limitations. 

We now specify the frame and overall geometry in which the final results will be inter- 
preted. We define the z' axis as the t direction in the ti rest frame. In the individual top 
quark rest frames, the W boson 3-momentum vectors (or equivalently the b quark vectors) 
can be denned relative to this axis by 

(sin #1 cos (j>i , sin 6\ sin , cos 9i ) , ( 1 ) 

(sin 62 cos <p2, sin 62 sin 2 , cos 9 2 ), 

where here W\ = W + and W2 = W~ . The angles 0i and 02 are invariant under boosts 
along the z' direction. In what follows we will be interested in the relative azimuthal angle 
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Wi 
W2 
W 2 



of the W bosons / b quarks in the ti frame about the z' axis and therefore define the angle 
$ = (02 — <Pi)- This angle represents the azimuthal separation of the two decay planes as 
indicated in Fig. |2J We note that although we have indicated that the W + /W~ bosons and 
the b/b quarks are distinguishable from one another, this could in principle not be the case 
experimentally. However, the results that follow still hold when the proper identification 
of states is made. What is relevant is that for each event a W boson and b quark can be 
paired uniquely to reconstruct a single top quark. In this case, we can designate either 
reconstructed top quark as t\, define z' = t\, and associate W\ with the W boson resulting 
from the decay of t\. The relevant angles are then once again defined according to Eq. ([I]) 
and displayed in Fig. [2] with the necessary replacements. 

We calculate the full matrix element directly for the process gg —> ti — >■ bbW + W~ using 
the diagrams shown in Figs. []Ja)-QTc) assuming real, on-shell W bosons. In addition, we 
make use of the narrow width approximation for the top quark propagators and integrate 
over all angles except $ = (02 — 0i) to find the normalized, parton- level result for the 



distribution in $: 

1 dlTQCD / 1 



<5qcd d$ \2n 



(2) 



where f3 t = a/1 — 4m|/s, p w = (Mn//m 4 ) 2 , and the coefficient is defined as 



r m - ( V-fiA / A(33-31^)-(l-^)(33-2/3 t 2 )tanh- 1 (A) - 
QCMW V Pt )\ A(59 - 31/3, 2 ) _ 2(33 _ m + ffl tanh -i (A) J 

The overall partonic level cross section is given by 



) [2 (33 - 18/3 2 + (3f) tanh- 1 ^) - A (59 - 31/3 2 ) ] (4) 
xBR(t->bW) 2 , 

where BR(t -)■ bW) = BR(t -> bW + ) = BR(t -> bW~) is the top quark branching fraction. 

We plot the coefficient Cqcd(A) in Figs. [3](a) and[3]^b) as a function of both f3 t and vl. 
In the latter, we take a top mass of m t = 174.3 GeV, which we use throughout the paper. 
As is evident from this figure, it is a well-behaved function which takes the limits 



°qcd (sj 



48 
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FIG. 3: Amplitude factor Cqcd as a function of (a) fit and (b) \f§. 



a 



QCD 



(5) 



1 as f} t -»• 
as fit ->• 1 

The cos $ dependence of the distribution is therefore maximal right above the production 
threshold of the ti pair and then rapidly goes away at higher energies. This behavior stems 
from the different energy dependences of the relative contributions to the total amplitude 
coming from initial state gluons of varying helicity, as discussed in Ref. {26 ]. 

We plot the resulting overall behavior of Eq. fl2]) in Fig. [4] for several different values of 
partonic center-of-mass energy vl- The strong dependence on y/§ is clearly demonstrated, 
with the distribution appearing nearly flat by V§ ~ 800 GeV. What we are now interested 
in is how contributions from New Physics may introduce deviations into this behavior. 



III. NEW PHYSICS MODEL 

As a simple and illustrative example, we consider a scalar particle 2 h coupling to the top 
quark according to the following Lagrangian: 

2 While we suggestively label the state as h, it is not necessarily associated with the Higgs boson of the 
Standard Model, as we intend our results to be more general and model independent. This choice is 
simply motivated by the fact that much of the formalism has been borrowed from previous studies more 
directly applicable to a light, SM-like Higgs boson. 
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if D - (— ) ht(AP L + A*P /? )t 



(6) 



where 



A = (^ + iA,) = |4e te - 



(7) 



We restrict a to the range — | < a < | without loss of generality. This gives a completely 
scalar (pseudoscalar) coupling when a = (±| ). In particular, when \A\ = 1 and a = 0, this 
is the same as the top quark coupling to the SM Higgs boson. Other values of a, however, 
would indicate CP violation. Such couplings can arise, for example, from mult i- Hig gs models 
with complex couplings after rewriting; the Lagrangian in the mass basis 27H30I] and have 



been studied in a variety of contexts 



12 



23 



3l| . Here we are primarily interested in probing 



the nature of such couplings using the angular correlations described above. 

We assume this scalar is heavy enough to decay on shell to a pair of top quarks. The 
existence of this scalar, then, introduces relevant new diagrams to the process gg — > tt — >■ 
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(a) 



(b) 



FIG. 5: Diagrams contributing to the process gg — > h — > tt — > bbW + W . Any colored fermions 
which couple to the new scalar boson h may appear in the fermion loop. 

bbW + W~ . These are given in Figs. |5]^a) and Et^b). These diagrams can interfere quite 



strongly with those of the purely QCD process, a point noted 



jy several authors when 



considering the effects of heavy scalars on gg — > tt production [15|, [18| . We will take special 
care to account for this interference analytically when discussing the complete azimuthal 
correlations below. 

Because of the scalar nature of h, the production process gg — > h^*' can be calculated 
separately and summed between the two contributing diagrams to give the amplitude: 



iM(gg -)• fc w ) = (i5 ab ) 



where tq = (Atjiq/s) and 



(£) (^X>(-«>) (-^fe 



P2) + P2P1) 



(8) 



Q 



Mr) 




Mr) 


= r/(r), 


f(r) 


-{-. 



(l + (l-r)/(r)) 



(9) 



arcsm 
In 



1 

HVi— F 



1 2 



VK 



if r > 1 



if r < 1 



l-Vl-T, 

The functions Ah(t) and Aa{t) are normalized such that Ajj{t), A^(t) — > 1 in the limit 
r — > 00 and are well known from searches of CP-even and CP-odd Higgs bosons produced 



via gluon fusion (see, for example, Ref. 32[ and references therein). For convenience we will 
now define 
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F H = ^Ah(tq), (10) 
Q 

F A = J2 A a(tq)- 

Q 

The sums appearing here and in Eq. ([8]) run over all the quarks that may contribute to the 
loop. It is assumed that any other quarks coupling to the scalar would do so according to 
Eq. ([6]) with the replacement m t — > vtlq. This also makes the assumption that the factor 
A is universal for all heavy fermions. While not necessarily generic, this will simplify the 
discussion below. We also note that in the above calculation, we have assumed that h could 
in principle be off shell such that it is not necessarily true that s = m^. Therefore, we use 
Tq = {Attlq/s) rather than the more commonly used value of tq = (Am^/rnfy seen in the 
literature. 

For completeness, we note that in the large r limit, the effective coupling of the scalar to 
two gluons can be characterized by the effective Lagrangian term 

=2f D --gugg h G^G^ - -gAg g h G^G^ , (11) 
where G a = |e M „ po -G p<T ' a is the dual to the gluon field strength tensor and 



v» = -(p A « F »' < 12 > 
3A » = " a,Fa ■ 

We note here that we use the definition e 0123 = (—£0123) = +1, which implies that the 
relation 



trh^VTV] = -4ze^ ff (13) 

was used in the calculation of gAgg- The overall and relative signs of the two terms in Eq. (JXXJ) 
[when substituting in the values from Eq. ( 1121) ] are often unimportant to the studies of Higgs 
searches. However, we will be concerned with the interference with the purely QCD process 
and therefore the proper signs are critical. Furthermore, when dealing with the case of CP 
violation - in which neither Ar nor Aj are zero - both terms will make a contribution to 
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this interference. For this reason we emphasize our use of e 0123 = +1, as a different choice 
would change the sign of gA gg in Eq. ffT2|) and of the CP-odd term in Eq. ([8]). 

We will consider two scenarios within this overall framework. In the first, we consider 
the top quark to be the only fermion coupling to the new scalar. The top quark is then the 
only colored fermion appearing in the loops in Fig. Eland in the sums in Eqs. ([H]) and (TTU|) . 
We refer to this as the "top-only" scenario. 

The other scenario we consider is one motivated by the possible existence of a fourth 
generation of fermions. There are then fourth-generation quarks that could also run in the 
fermion loops along with the top quark. In addition to the boost in the production cross 



section these fermions would provide [321 ]. such a scenario could also relieve the tension a 
heavy scalar may have with electroweak precision observables if it were to be associated with 
the Higgs boson 



33 



34J. For our purposes, we will only assume for this scenario the existence 
of two additional quarks with mq = 500 GeV but will, for simplicity, refer to this as the 
"fourth-generation" scenario. This particular value of mq is chosen not only for consistency 
with direct search bounds on fourt h-g eneration quarks - with currently published results 



placing a limit of my > 350 GeV [35l . l36| and preliminary results indicating an exclusion 
m t' ^ 450 GeV p\- but also to illustrate the effects of contributions to Eq. ([TO]) with tq > 1 
and therefore which provide contrasting behavior to the top-only scenario, where Tj < 1. 



We note that values of mq much 



unitarity of the scalar sector 



38 



arger than this could run into problems with perturbative 



39|. 



In both scenarios, we assume that h — > ti is the only decay mode with an appreciable 
width, i.e., the partial widths h — >■ gg and h — > 77 are ignored and any new colored fermions 
are assumed to have mq > m^/2. The total scalar decay width is then given by 

Unless otherwise stated, we fix the magnitude of A from Eq. ([7]) to be \A\ = 1. We note 
that the observability of the scalar signal depends very strongly on the above assumption. 
If, for instance, the scalar is treated as a very heavy Standard Model Higgs boson with no 
enhancement of the cross section, the implied increase in the scalar decay width from direct 



couplings to weak gauge bosons would wash out any observable effects [2l| . The width above 
represents the smallest possible value for a scalar coupling to the top according to Eq. ([6]) 
and therefore can be thought of as a limiting scenario. 
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IV. TOP PAIR INVARIANT MASS DISTRIBUTIONS 



Prior to our discussion of the effects such models have on the azimuthal correlations of 
the decay planes of the top quarks, we will first discuss the relative contributions to the 
invariant mass of the top pair coming from QCD, the scalar signal, and the interference 
between them. The primary reason for this is that the deviations from the QCD description 
of the azimuthal correlations can only in principle be observed in small regions of the top 
pair invariant mass where the signal is localized. Furthermore, the size of the deviations 
will depend on the relative strength of the signal to the QCD background and therefore it 
is important to first have a thorough understanding of this distribution before proceeding. 

At LO, the parton-level cross section for gg — >■ h — >• ti —> bbW + W~ and its interference 
with QCD are given by 




x BR{t->bW) 2 , 



where we define 



C H = Re(F H 0» h ) 
C A = Re(F A &' h ) 



(17) 



with Fjj and F A from Eq. ([TO]) and denote the Higgs propagator by 



(a - m 2 h ) + im h Y h 



(18) 



These are calculated in the narrow width approximation for the top quark propagators 



and are in agreement with previous calculations in the CP-conserving limits 



15|. The 
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term ^'h, which appears in Eq. (fT7|) . is not the full scalar propagator, but the scalar 
propagator stripped of an overall factor of i. It is important to note that unlike in the 
signal cross section in Eq. ( II 5p . this factor appears in the interference term in Eq. (FTO]) 
unsquared via the functions Ch and C\. This results in cxint(s) switching signs near, but 
not typically at, s = m^. This is an important effect in the overall results, as the magnitude 
of the interference term can be rather large relative to that of the signal. The resulting 
negative contributions can give rise to a "peak-dip" structure in the top pair invariant mass 



distribution, whereby the final resu 



with comparable magnitude 







15 



t may lie both above and below the QCD prediction 



18| 



The total parton-level cross section for gg — > tt — >■ bbW + W is the sum 

Otot = OqCD + (Th + 0"INT • (19) 

Each of the formulas in Eqs. @, f|T5|) . ffT6]) . and f lT9|) gives the results for the process gg —> tt 
(absent the decay of the top quarks) with the removal of the BR(t — > bW) factors. We can 
then integrate over a set of parton distribution functions (PDFs) to find the invariant mass 
distribution of the top pair at the LHC using 



J dyf g / P (Vre y ,n F ) f g/p (^e y , fx F ) a(r) , (20) 



da 

where f g / p is the gluon PDF, fi F is the factorization scale, y is the rapidity of the center-of- 
mass system in the lab frame, and r = (s/s). Then 



da da 



dM ti dV§ 
The total cross sections are then given by 




(s/s)) . (21) 



a = J dydrf g/p (y/re y , /j f ) f g/p [yfre y , fj, F ) a(r) . (22) 

We plot the invariant mass distribution of the top pair from gg — > tt in Figs. |6]^a)-[6]T) 
for QCD and for the top-only and fourth-generation scenarios we consider. We use several 
scalar masses ranging from m h = 400 GeV to 900 GeV and the phases (a/vr) = 0,0.25, 
and 0.50. In calculating these distributions, we have used a total center-of-mass energy 
-y/i = 14 TeV, CTEQ6L1 PDFs for proton-proton collisions 40J, and fixed renormalization 
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FIG. 6: Invariant mass distributions of the top pair from the process gg — > tt at the LHC with 
y^s = 14 TeV for the top-only (a)-(c) and fourth-generation (d)-(f) scenarios. These are calculated 
for several different values of scalar masses and for the phase a. The magnitude of the normalized 
scalar-top coupling A is fixed at \A\ = 1. 
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and factorization scales of /ir = [ip = 350 GeV~ 2m t . In addition, we have used one- 
loop running of a s evaluated at /i^. These results are in good agreement with those of 



Refs. 



15 



17 



20]. 



When comparing the two scenarios for a fixed value \A\ = 1, it is clear that for a given 
scalar mass and phase a, the fourth-generation scenario gives a much larger deviation from 
the QCD prediction. This is due to the well-known enhancement to the scalar production via 
gluon fusion resulting from the existence of new heavy quarks running in the quark loops 
of Fig. [5] 32|. Indeed, while the fourth-generation scenario provides a scalar signal that 
appears observable in principle up to the largest mass considered here (rrih = 900 GeV), the 
top-only case does not have easily observable signals much above the top quark production 
threshold to around rrih « 500 GeV. This would not be much improved for larger values of 
\A\, as this tends to simply broaden the distributions seen in Figs. |6^a)-|6]T). 

For a given scenario and given scalar mass, the effect of a variation in a on the top 
pair invariant mass distribution is twofold: as seen in Eqs. ( TT5l) and ( JT6l) . there is a larger 
threshold suppression of both the signal and interference contributions to the total cross 
section asa->0 (the case of a CP-even scalar), implying smaller overall signals relative to 
the case for when a — > 0.5n. However, there is a corresponding suppression in the scalar 
decay width of Eq. (|T4|) . For smaller scalar masses, the result is larger, narrower peaks in the 
invariant mass distribution as a — > 0, seen most dramatically in the case of rrih = 400 GeV. 
At higher scalar masses, however, larger peaks are seen in the limit a — > ±0.57r. 

As an example showing the direct relationship between the different contributions to the 
invariant mass distribution, we plot the signal, interference, and QCD terms for a scalar 
mass rrih = 800 GeV and a = 0.257T in the fourth-generation scenario in Fig. [71 As discussed 
above, the interference term is non-negligible when compared to the signal and switches sign 
for a value of \/I = M t i that is somewhat less than the location of the signal resonance at 
y/§ = rrih. This results in a peak of the total distribution that is shifted from the resonance, 
followed by a modest, extended dip relative to the purely QCD case. The negative values 
taken by the interference term are important to note, as they serve not only to reduce the 
total cross section in these regions but will also result in a phase shift of it in the interference 
contribution to the normalized azimuthal distributions, a point we will return to below. 
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FIG. 7: Invariant mass distributions of the top pair from the process gg — > ti at the LHC with 
y/s = 14 TeV for the fourth-generation scenarios with a scalar mass rrih = 800 GeV and phase 
a = 0.257T, showing the individual signal, QCD, and interference contributions to the overall 
distributions. The magnitude of the normalized scalar-top coupling A is fixed at \A\ = 1. 

V. AZIMUTHAL CORRELATIONS 

Following the same procedure described in Sec. [Ill and making use of the same event 
geometry, we calculate the normalized, partonic distributions for the angle $ = {<p2 — <j>i) 
for the signal and interference terms separately. These are given by 



1 dan 



1 

2^ 



1- - 



tt\2 (1 - 2p w f {A 2 R f3 2 t - A]) cos$ + (2AjA R (3 t ) sin$ 



47 (l + 2p, 



+ A}) 



(23) 



and 



1 da 



INT 



(Tint d& 



1 

2^ 



7T\2 

4. 



2p w f \ (2A\$C H - 3A]C A ) 
2 Pw f \ (2Al(3fC H + 3A 2 C A ) 



cos $ 



(24) 



(A I A R p t )(3C A + 2C H ) . 
(2A R ffC H + 3AjC A ) Sm 
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where Cha and defined above. These can be greatly simplified and more easily 

interpreted if we now define phases x an d £ according to 



.4, 



{A R p t + iAj) 



\A x \e lx 



These are related to the phase a, the main parameter of interest, by 



tanx 
tan£ 



Aj 

PtA R " 
SAjC A 



— I tana 

Pt, 



3 Re {Fa&' \, 



2A R (3 t C H \2Re(F H 0» h ] 
With the above definitions, Eq. (l23~l) and (EH become 



— I tan a 

Pt, 



1 d&u 



2n 



4/ 1 



2/^ 



2p, 



cos($ - 2%) 



and 



(25) 



(26) 



(27) 



1 da 



INT 



(Tint ^<3> 



1 

2^ 



As discussed in Refs. 



22 



Try (l-2p w ) 2 cos($-(x + Q) 
4 V (l + 2p^) 2 cos( X -0 



(28) 



4l|, the phase x appearing in Eq. fl27j) gives a rather direct 



experimental probe of the CP nature of the scalar particle and the phase a when considering 
the signal distribution in isolation: for highly boosted tops, and therefore large values of /3t, 
and there is a phase shift in the $ distribution of ~ —2a. The precise behavior of x 
as a function of y/§ is given for several different values of a in Fig. [HJ We note that in the 
CP-even (a = 0) and CP-odd (a = 0.57r) cases, x takes the same value as a independent of 
f3 t . In the CP-odd case, there is then a phase shift of ir and the signal distribution is exactly 
in phase with the QCD-only distribution, while in the CP-even case, the two distributions 
are correspondingly out of phase by 7T. 

What is more interesting, however, is the possibility to detect CP violation by means of 
a small phase shift in this distribution, a quantity which is often measured using tools such 
as "up-down" (or similar) asymmetries Q, [l8|, essentially comparing the values for $ > 
to those with $ < 0. While the full $ distribution is perhaps more difficult to measure than 
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FIG. 8: The value of the phase x as a function of ys for several different values of the phase a. 
The dotted lines give the corresponding fixed values of a. For a < 0, the resulting behavior is the 
mirror image about the \/I axis. 

an asymmetry - if only for purely statistical reasons - it is interesting in that it contains 
more information: not only does it distinguish between the CP-even and CP-odd cases, 
which an up-down asymmetry does not, but it shows how simply and directly the violation 
is related to the phase a. 

It is also interesting to note that the normalized signal distribution at the parton level 
is completely characterized by the one parameter a through the phase x\ it carries no 
information of the magnitude \A\, the total width Th, nor of any other details concerning 
the production of the scalar itself. This is due to the scalar nature of the signal: because 
there are no spin correlations in the production and decay of h, the normalized $ distribution 
for cr(gg — > h — > tt — > bbW + W~) is the same as that for the partial width T(h — y ti — > 
bbW + W~). Normalized thus, Eq. fl27|) is a fairly generic result for scalars which couple and 
decay to weakly decaying particles via a coupling like Eq. fl6]), a point noted in (41]. 

While the normalized distribution of the signal can be completely characterized by the 
phase the distribution derived from its interference with the QCD process also depends 
on the phase £. Like Xi the phase £ takes on the same value as a for purely CP-even or 
CP-odd couplings. In these cases - and more generally for cases in which £ = x ~ the 
normalized distribution derived from the interference term takes the exact same form as 
that for the signal and will therefore serve only to enhance or minimize the effect of the 
signal depending on the sign of cxint- However, in the CP-violating case the behavior of £, 
and therefore the effect of the interference on the $ distribution, is more complicated. We 
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FIG. 9: The value of the phase £ as a function of \/§ for the choice a = 0.257T and several values 
of the scalar mass for the (a) top-only scenario and the (b) fourth-generation scenario. This is 
compared to the value of x-, a ^ so as a function of plotted on the same scale and the fixed value 
of a = 0.25vr. 

plot the value of £ as a function of vl in Figs. |9^a) and|9^b) for the fixed value a = 0.25-7T 
and for several different values of the scalar mass m^. This is done for both the top-only 
and the fourth-generation scenarios. In the context of its use in Eq. ff28l) . we can restrict 
£ to the range —n/2 < £ < n/2 without loss of generality. This will be assumed for the 
remainder of this paper. 

For a given scalar mass, there is a region in both the top-only and fourth-generation cases 
where £ switches sign, taking on its full range of values in the process. This occurs near, but 
not generally at, a/I = mh and is again related to the pole in the propagator appearing in the 
functions Ch and Ca- The complications that would arise from this are generally minimized 
by the fact that the interference contribution to the partonic cross sections is switching sign 
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here as well and is therefore in the region where it contributes least to the overall behavior 
of the full distribution. What is more important, then, is the behavior of £ outside this 
region, where the interference plays a much larger role. It is here that the two cases (top- 
only and fourth-generation) are quite different. Outside of the pole region, the value of £ 
in the fourth-generation case takes a relatively stable value lying slightly apart from \ for 
all scalar masses, while no such behavior is seen for the top-only case. This is rooted in 
the different contributions that the top quark and the fourth-generation quarks make to the 
functions Fh and Fa- The new heavy quarks have tq > 1 such that their contributions are 
purely real while the top has r t < 1 and therefore gives a complex contribution. The former 
contributions dominate in the fourth-generation case, minimizing the effect of the complex 
part of the scalar propagator, and therefore the scalar width, in Ch,a- This results in the 
ratio Ca/Ch ~ 1 in Eq. ff26l) for lower values of which then grows slightly approaching 
the tq = 1 threshold at Vs = 1 TeV. The main conclusion that can be drawn here is that 
in the fourth-generation case, the interference term will give contributions to the overall $ 
distribution that are approximately in phase with the pure signal contribution as long as 
<3"int(s) > 0, while in the top-only case this is much more dependent on the scalar mass and 
the particular top pair invariant mass region of interest. 

Now, given all the individual contributions to the $ distribution, the complete observable 
behavior will be given by 



where each at signifies that the parton-level results have been integrated with PDFs over 
some region of top-pair invariant mass. We note here that for each contribution we can 
define the normalized distribution and relative weight as 



1 da tot 
otot d§ 




(29) 



/<(*) 



1 doi 



(30) 




(31) 



such that the total distribution can be given as 



/«($) = (w QC d/qcd($) +%/hW +wint/int($)) , 



(32) 
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with Yli w i = The total <J> distribution is therefore the sum of the individual normalized 
contributions weighted by their contribution to the cross section 3 . Like the individual terms, 
we can write 

/«(•) = (^) (l + C M Q 2 l^cosf* + W ) , (33) 

where C to t and 5 fo4 depend on the values of the various parton-level angular distributions 
integrated with the PDFs. We note, however, that when integrated over a small enough 
range of vl = M t i, these distributions should be a good estimate of the resulting behavior, 
as the PDF integrations cancel out if the weights and angular distributions are assumed 
to have little variation in that region. We can then make the statement that what we are 
ultimately interested in, then, is if deviations can be seen in one or both of Ctot = Cqcd 
and S to t = due to the effect of a scalar particle coupling to the top quark as described. 

In Figs. [T07 a)-[T07f). we plot the total normalized $ distributions in addition to the 
contributions from QCD, the scalar signal, and their interference at the LHC at y/s = 14 TeV 
for a scalar mass of rrih = 400 GeV and for (oc/tt) = 0,0.25,0.50. These distributions have 
been integrated with PDFs over the range 390 GeV < < 400 GeV. This particular range 
is motivated by the desire to avoid regions above the resonance in which the interference 
gives a negative contribution, thereby weakening the overall strength of the deviation. These 
plots are meant to give an example of the maximum effect the scalar can have on the <3> 
distribution, as = 400 GeV gives some of the largest signal cross sections possible, as seen 
in Fig. E^a)-|(3T). Rather large departures from the QCD prediction are seen, in particular, 
for the fourth-generation scenario. For instance, in the CP-even case (a = 0), the resulting 
distribution is completely out of phase with the QCD prediction, albeit with a reduced 
amplitude. This scenario is rather unrealistic when considering this particular scalar mass, 
however. Such a large spike in the invariant mass distribution of top pairs would likely have 
been observed in the early running of the LHC, if not at the Fermilab Tevatron, as this is 
precisely the region with the largest statistics. 

We show a more realistic possibility in Figs. HTf a)- IITlf f). This time we use a scalar mass 

3 We note one slight subtlety here: it can be shown that while u>int —> as (\ — £) —> ±7r/2, there will 
generally be a small non-zero contribution to ftot(&) from the interference term here due to the factor 
cos _1 (x — £) appearing in Eq. P8|). 
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FIG. 10: Total normalized $ distributions for the process gg — >■ ti — )• bbW + W~ , in addition to 
their individual contributions from QCD, the scalar signal, and their interference, at the LHC with 
y/s = 14 TeV for the top-only (a)-(c) and fourth-generation (d)-(f) scenarios. These are calculated 
for the scalar mass = 400 GeV and for the phase (a /it) = 0, 0.25, 0.50 and have been integrated 
with PDFs in the range 390 GeV < \f§ < 400 GeV. The magnitude of the normalized scalar-top 
coupling A is fixed at \A\ = 1. 
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FIG. 11: Total normalized $ distributions for the process gg — >■ ti — )• bbW + W~ , in addition to 
their individual contributions from QCD, the scalar signal, and their interference, at the LHC with 
yfs = 14 TeV for the top-only (a)-(c) and fourth-generation (d)-(f) scenarios. These are calculated 
for the scalar mass = 800 GeV and for the phase {ct/n) = 0, 0.25, 0.50 and have been integrated 
with PDFs in the range 790 GeV < \/I < 800 GeV. The magnitude of the normalized scalar-top 
coupling A is fixed at |^4| = 1. 
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FIG. 12: The data points are a "binned" version of the full expected $ distribution in the region 
790 GeV< V§ <800 GeV for a scalar mass = 800 GeV and phase a = 0.257T at the LHC with 
yfs = 14 TeV. The error bars signify the expected magnitude of the statistical uncertainty in each 
bin for a luminosity of (a) 100 fb -1 and (b) 500 fb" 1 . The dashed red curve gives the LO QCD 
background prediction for distribution from the process gg —> ti —> bbW + W~ . 

of nih = 800 GeV and integrate over the region 790 GeV < \/i < 800 GeV for the same 
reasons described above. Not only are the signals less pronounced here in the invariant 
mass distribution, but the QCD contribution is roughly an order of magnitude smaller 
than in the previous case as well. The distributions in Figs. [TTT a)- fTTT f ) are therefore more 
representative of possible signals that may currently be hidden. The top-only case creates 
essentially no deviation from the QCD prediction, as expected from its lack of observable 
signal in the M t i distribution. However, in the fourth-generation case the deviations are 
small but seemingly distinguishable. This is aided by the fact that at this center-of-mass 
energy, the $ distribution predicted by QCD is nearly flat. Also due to this fact, the phase 
of the total distribution very nearly matches that of the signal, with 5 to t ~ —^X ~ — 2a, and 
the distributions for each of the three values of a shown are clearly distinguishable. 

It is important to note that the size of the deviations described above are only meaningful 
when understood within the context of the relevant systematic and statistical uncertainties 
of the experiment. While we do not perform a full Monte Carlo simulation, we attempt 
to get an estimate of the statistical uncertainties involved by plotting a "binned" version 
of the total normalized $ distribution. This is performed by integrating over regions of 
A($/7r) = 0.2, dividing by the bin size, and plotting the error bars that would be associated 
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with the counting statistics on the number of events in each bin predicted for a given 
luminosity. We do this in Figs. IT27 a) and H~27 b) for integrated luminosities of Jzf = 100 fb" 1 
and «if = 500 fb _1 , respectively, for the fourth-generation scenario with = 800 GeV 
and a = 0.257T. The QCD prediction is plotted as a smooth curve. What can be seen in 
the comparison between Figs. [T2T a) and [T2T b) is that in this more realistic scenario, a very 
large amount of data (on the order of hundreds of fb" 1 ) will be necessary to begin to have 
sensitivity to the deviations at this level, but that they are in principle possible to observe. 
For smaller scalar masses this requirement would be alleviated somewhat due to the more 
favorable statistics coming from more pronounced signals and larger overall rates. 

We close with a few comments on additional theoretical concerns. Firstly, it has been 



noted previously in Refs. 
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26| that the azimuthal angular correlations of the IV-boson 



decay products exhibit a stronger modulating signal than those of the W bosons themselves 
due to a 100% spin correlation of the down-type decay objects with the polarization of the 
parent top quark. While we have worked in the on-shell W approximation primarily for 
ease of calculation, we expect little improvement from such an extension of the calculation, 
as both the signal and background will be affected and therefore the question of detecting 
a deviation remains. A larger overall modulation could, however, improve the statistics 
necessary for detection and is therefore worth consideration. 

Secondly, we note here that we have dealt only with the gluon initiated process and 
have ignored the processes qq — > ti. This should be a fairly good approximation: while 
qq — > ti makes up a small but non-negligible fraction to the total pp — > ti cross section 
at -y/i = 14 TeV, it provides a continuum distribution much like the gg —> ti contribution 
and should therefore be insignificant when considering the small, isolated ranges of M t i used 
here. 

Finally, it is well known that higher-order contributions to the total pp — > ti cross section 
in QCD are quite large, resulting in K corrections of ~ 1.5 2lJ. These, of course, would have 
an effect on the da/dM t i distribution of both the QCD process and the signal. Any difference 
in the corrections between the two would change the relative weights used in the total angular 
distribution and therefore have an effect on the final results. In addition, it is possible that 
there could be modifications to the normalized angular distributions themselves. While all 
of these effects would be necessary for a more accurate understanding of any deviations 
that may be observed in these distributions, we believe the results here will serve as a good 
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approximation: because most higher-order corrections are related to initial- and final-state 
radiation, and because these states are the same for the QCD contributions and the scalar 
signals described here, we would expect them to be of comparable magnitude and therefore 
approximately cancel out in the normalized distributions we consider. In addition, the next- 
to-leading-order corrections have been calculated for the azimuthal angular distribution of 
the final-state leptons resulting from the W^-boson decays in the QCD process and these 



have been shown to be small at the LHC 



42]. 



VI. CONCLUSIONS 



We have presented the LO prediction for the distribution of the azimuthal angle between 
the decay planes of the top quarks in the process gg — >■ ti —> bbW + W~ . We have then 
sought to demonstrate the deviations that may arise in this distribution from possible New 
Physics contributions by specifically considering the case of a heavy scalar particle with a 
coupling to the top quark such that it can contribute via an on-shell decay. We investigated 
this possibility within the context of two separate scenarios: one in which this coupling to 
the top is the only relevant coupling of the scalar, and another, inspired by fourth-generation 
models, that includes couplings to two new, heavy quarks. We have placed special emphasis 
on understanding the contribution to the total azimuthal distribution from the interference 
of the signal with QCD and on the possibility of CP-violating couplings of the scalar to the 
top quark. We have shown that, in certain scenarios, large amounts of integrated luminosity 
may allow this distribution to be used to directly probe a, the phase associated with these 
couplings. 
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